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JLQCD collaboration recently started the Nf = 3 QCD simulations with the 0(a)-improved Wilson fermion 
action employing an exact fermion algorithm developed for odd number of quark flavors. It is found that this 
theory has an unexpected non-trivial phase structure in the (/3, k) plane even at zero temperature. A detailed study 
is made to understand the nature of the observed phase transitions and to find the way of avoiding untolerably 
large lattice artifacts associated with the phase transition. 



1. Introduction 

Including the dynamical quark loop effects in 
large-scale QCD simulations is one of the most 
important and urgent problem. While in the 
real world there are three light quarks, most of 
recent studies have concentrated on the case of 
two flavors of Wilson- type fermions because these 
systems can be simulated by the exact Hybrid 
Monte Carlo (HMC) algorithm. Recently, how- 
ever, a major progress has been made in the ex- 
act fermion algorithm applicable to odd number 
of flavors as reviewed in Ref. B. In particular, 
we have succeeded in developing an efficient algo- 
rithm for Nf = 3 QCD with the 0(a)-improved 
Wilson fermion action in the framework of the 
HMC algorithm. Leaving algorithmic details to a 
separate report || , we present here the first result 
toward realistic Nf = 3 QCD simulations with 
the 0(a)-improved Wilson fermion. We found 
that this theory has an unexpected first-order 
phase transition at zero-temperature, which gives 
a strong constraint on the form of lattice actions 
suitable for large-scale simulations. 



2. Phase structure for plaquette gauge ac- 
tion 

The partition function we study is defined by 
Z= [vU(det[D ud }f(det[D s ])e- s ^ u \ (1) 



Here S g (U) is the gluon action 



s g (u) = f [■ 



c y^ lUixi +ci 
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where Co = 1 — 8ci and Wi x i and Wi X 2 are the 
plaquette and rectangular Wilson loops, respec- 
tively, and (det[Z? U( 2]) 2 represents the contribu- 
tion from degenerate u and d quarks whereas s 
quark effect is given by (det[£) s ]) with D q (q = 
ud, s) the O(o)-improved Wilson-Dirac operator. 

We start our analysis employing the plaque- 
tte gauge action (ci = 0), determining the clover 
coefficient c sw by tadpole-improved one-loop per- 
turbation theory; 
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* Presented by M. Okawa. 



with (P) evaluated in the quenched approxima- 
tion at the corresponding value of (3. 

Given the lattice action, we have three tunable 
parameters, the gauge coupling /?, the hopping 
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Figure 1. Thermal cycles on an 8 3 x 16 lattice at 
(3 = 4.6, 4.8, 4.85, 4.9, 4.95, 5.0, 5.05, 5.1, 5.15, 
5.2, 5.25, 5.3, 5.4, 5.6, 5.8, 6.0 from bottom to 
top. 

parameter K u d for ud quarks and k s for s quark. 
For simplicity we set k uc i = k s = n in this work. 
To see the global phase structure in the (/?, k) 
plane, we perform rapid thermal cycles on 4 3 x 
16 and 8 3 x 16 lattices. Fixing (3 we increase k 
from k = in steps of 0.01 until the algorithm 
fails to converge and then decrease k to 0. The 
procedure is repeated for several values of (3. In 
this analysis, we omit the noisy Metropolis test 
for the correction factor de^H^o] ||] and fix the 
order of the Chebyshev polynomial to n = 100. 

The result of the thermal cycle on an 8 3 x 16 
lattice is shown in Fig. [|. We clearly observe hys- 
teresis loops at (3 = 4.95 and 5.0. On a 4 3 x 16 
lattice, the hysteresis loop is seen in a wider range 
of j3 at 4.8 < (3 < 5.1. To understand the nature 
of the hysteresis in more detail, we perform sim- 
ulations starting from both ordered and disorder 
configurations with fixed values of (3 and k on 
8 3 x 16 and 12 3 x 32 lattices. Some examples are 
shown in Fig. ^. Here we make simulations exact 
by including the correction factor det [Woo]. We 
observe clear two state signals for several sets of 
parameters, which demonstrates that the hystere- 
sis in the thermal cycles is caused by a first-order 
phase transition. 

In Fig. H we plot the location of the observed 
phase transition line (n = nx((3)) on the (/?, k) 
plane for 8 3 x 16(filled circles) and 12 3 x 32 (open 
upper-triangles) lattices. Also shown in Fig. || are 
the location of the hysteresis observed on 4 3 x 16 
and 8 3 x 16 lattices (open squares and circles). 



0.52 
0.50 
0.48 
0.46 
0.44 
0.42 




>c=0.1350 



>c=0.1345 

/ \ / 



500 



1000 
Traj. 



1500 



2000 



Figure 2. Two-state signals on a 12 3 x 32 lattice 
at /3=4.88 and c sw =2.15. 



0.14 



* 0.13 



0.12 



w 



v- 



▼ w 



Q 



-B- 4 x16,k- x 

>&* 8 3 x16,k x 

+ 8 3 x16,k x , exact 

>*< 12 3 x32,k x , exact 

^ 8 \ 16.x; 

4.6 4.8 

j8 



5.0 



5.2 



Figure 3. Phase diagram in ([3, n) plane. 

The position of the phase transition line signifi- 
cantly moves from 4 3 x 16 to 8 3 x 16. However, 
it stays in the same place when the lattice size 
increases from 8 3 x 16 to 12 3 x 32. In Fig. ||, we 
plot the plaquette values in the two phases at the 
phase transition line as functions of (3 for 8 3 x 16 
and 12 3 x 32 lattices. The gap in the plaquette is 
almost independent of the lattice size. These find- 
ings strongly suggests that the first-order transi- 
tion line persists in infinite lattice volume and 
zero temperature. 

Figure || also shows that the gap in the plaque- 
tte decreases almost linearly in (3 and vanishes 
at around (3 — 5.0. Since no signs of transi- 
tions are observed at (3 > 5.0 in Fig. [l], we do 
not think that this point corresponds to a tricrit- 
ical point where a first-order phase transition line 
changes to a second-order one which then extends 
toward weaker couplings. Rather, the observed 
first-order transition must be a lattice artifact re- 
stricted to strong coupling regions. 
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Figure 4. Plaquette values in two phases at the 
first-order transition line as functions of [3. 
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Figure 5. Thermal cycle with RG gauge action at 
,3=1.50-2.25 in steps of 0.05 (bottom to top). 

3. Phase structure for improved gauge ac- 
tion 

To confirm this possibility, we repeat thermal 
cycle simulations employing the improved gauge 
action with c\ = —0.331 (RG action ||) and 
ci = — 1/(12(P)5) (tadpole-improved Symanzik 
action []4|). Here c sw is again determined by 
tadpole-improved one-loop perturbation theory. 
Figure || shows the result of thermal cycles for the 
RG action on an 8 3 x 16 lattice. No signs of hys- 
teresis loops are seen, and hence non-trivial phase 
structure is absent with this gauge action. Hys- 
teresis is also not seen for the tadpole-improved 
Symanzik gauge action. We conclude that the 
first-order phase transition observed for the pla- 
quette gauge action has the nature of a lattice 
artifact. 

4. Implications 

It has been known [|| that pure SU (3) gauge 
theory with fundamental and adjoint couplings 



has a first-order transition for positive values of 
the adjoint coupling. It is an interesting possi- 
bility that the clover term produces a positive 
adjoint coupling in an effective theory after inte- 
grating out the fermion degrees of freedom, which 
eventually causes the first-order phase transition 
observed in this work. In fact we have checked 
that there are no sign of phase transitions in 
thermal cycles made with the unimproved Wil- 
son fermion action (c sw = 0) and plaquette gauge 
action (ci = 0) on a 8 3 x 16 lattice, which sup- 
ports our interpretation. 

Our finding poses a serious practical problem 
whether realistic simulations with reasonable val- 
ues of the lattice spacing a are possible with the 
0(a)-improved Wilson fermion and the plaquette 
gauge action. To study this problem, we evalu- 
ate the lattice spacing from static quark poten- 
tial at (3 — 5.0. Our results are a -1 = 1.53(2) 
GeV (k = 0.1320), 2.06(5) GeV (n = 0.1330) and 
2.58(4) GeV (k = 0.1338). Continuum extrapola- 
tions using lattices with a -1 starting at 2.58 GeV 
is not practical. Our natural choice, then, is to 
use improved gauge actions for realistic Nf = 3 
QCD simulations with the (9(a)-improved Wilson 
fermion. 
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